AD-A140  980 


Mechanics  and  Materials  Center 
TEXAS  A&M  UNIVERSITY 
College  Station,  Texas 


THE  DYNAMIC  STEADY-STATE  PROPAGATION  OF  AN 
ANTI-PLANE  SHEAR  CRACK  IN  A  GENERAL  LINEARLY 
VISCOELASTIC  LAYER 


JAY  R.  WALTON 


s 

T 


DTIC 

ELECTEI 
MAY  1  1  1984 


•  *  **•/*•> 


MM  4867-84-4 


OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 

CONTRACT  NO.  N00014 -83 -K-0211 


FEBRUARY  1984 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


84  O' 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  This  PAGE  (Wh»n  Dmtt  Enfrnd) 

f  REPORT  DOCUMENTATION  PAGE 


11.  REPORT  NUMBER 


IIAC  READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


1 4.  TITLE  (and  Subtitle) 


'Hie  Dynamic  Steady-State  Propagation  of  an 
Anti-Plane  Shear  Crack  in  a  General  Linearly 
Viscoelastic  Layer 


7  AUTHORS 


Jay  R.  Walton 


5  TYPE  OF  REPORT  &  PERIOD  COVERED 

Technical  Report 


6.  PERFORMING  ORG.  REPORT  NUMBER 

VM  4867-84-4 


e.  CONTRACT  OR  GRANT  NUMBERfs) 

N00014-83-K-0211 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


».  PERFORMING  ORGANIZATION  NAME  ANO  AODRESS 

Texas  A&M  University 
Mechanics  &  Materials  Center 
Colleqe  Station.  TX  77843 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
Department  of  the  Navy 

Arlington.  VA  22217  _  _ 


4.  MONITORING  AGENCY  NAME  ft  ADORESSOt  different  from  Controlling  OZ/ire;  15.  SECURITY  CLASS,  (of  thm  report  j 


NR  064-520 


12.  REPORT  DATE 

February  1984 


13.  NUMBER  OF  PAGES 

15 


Unclassified 


IS*.  DECLASSIFICATION  DOWNGRADING 
SCHEOULE 


16.  DISTRIBUTION  STATEMENT  (of  thla  Report) 


DISTRIBUTION  STATEMENT  A 

v  Approved  lot  public  releoMl  s 
I  *  Distribution  Unlimited  -sm 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 


*9  KEY  WORDS  ( Continue  on  reverse  aide  if  necessary  and  identity  by  block  number) 

Dynamic  viscoelastic  crack  propagation 
Dynamic  stress  intensity  factor 
Riemann-Hilbert  boundary  value  problon 


ABSTRACT  .  Continue  nn  reveres  side  If  necessary  and  Identify  hv  block  number) 

Ain  a  previous  paper,  the  dynamic,  steady-state  propagation  of  an  semi- 
infinite  anti -plane  shear  crack  was  considered  for  an  infinite,  general 
linearly  viscoelastic  body.  Under  the  assumptions  that  the  shear  modulus 
is  a  positive,  non- increasing  continuous  and  convex  function  of  time, 
convenient,  closed-form  expressions  were  derived  for  the  stress  intensity 
factor  and  for  the  entire  stress  distribution  ahead  of  and  in  the  plane  of 
the  advancing  crack.  The  solution  was  shown  to  have  a  simple,  universal 


_ Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGF  ■»‘hrn  D»r«  VTTi 


DD  ,  1473 


SECURITY  CLASSIFICATION  OF  This  RA6Cf»h«n  Data  Batata4) 


>  dependence  upon*  the  shear  modulus  and  crack  speed  from  which  qualitative 
information  can  readily  be  gleaned.  Here,  the  corresponding  problem  for  a 
general,  linearly  viscoelastic  layer  is  solved.  An  infinite  series  repre¬ 
sentation  for  the  stress  intensity  factor  is  derived,  each  term  of  which 
can  be  calculated  recursively  in  closed-form.  As  before,  a  simple  universal 
dependence  upon  crack  speed  and  material  properties  is  exhibited. 


i 


S 


DTIC 

ELECTE| 
“AT  1  1 1984 


Unclassified 


2 


Problem  Formulation 

The  analysis  of  dynamic  fracture  models,  that  is  models  based  upon 
the  equations  of  motion  of  continuum  mechanics  rather  than  the  equations 
of  equilibrium,  has  received  considerable  attention  recently  in  the 
applied  mechanics  literature.  Principally  for  the  sake  of  mathematical 
convenience,  most  of  these  studies  have  been  in  the  context  of  linear 
elasticity.  Consequently,  for  elastic  material,  many  canonical  boundary 
value  problems  have  been  solved,  either  in  closed  form  or  numerically, 
for  both  steady-state  and  transient  modes  of  crack  propagation.  For 
viscoelastic  material,  however,  much  less  progress  has  been  made  in 
constructing  convenient  analytical  solutions  to  even  the  simplest 
dynamic  crack  problems. 

In  [2] ,  an  analysis  was  presented  of  the  dynamic,  steady-state 
propagation  of  a  semi-infinite,  anti-plane  shear  crack  in  a  general, 
infinite,  homogeneous  and  isotropic  linearly  viscoelastic  body.  (See 
[2]  for  a  discussion  of  other  relevant  studies  of  dynamic  viscoelastic 
crack  propagation.)  With  only  very  weak  assumptions  on  the  shear 
modulus  (specifically,  that  it  be  a  continuously  differentiable,  convex 
and  non- increasing  function  of  time),  simple  closed  form  expressions 
were  constructed  for  the  stress  intensity  factor  and  the  entire  stress 
field  in  front  of  and  in  the  plane  of  the  advancing  crack.  Moreover,  it 
was  shown  that  if  v,  c  and  c*  denote  the  crack  speed  and  elastic  shear 
wave  speeds  corresponding  to  the  zero  and  infinite  time  values  of  the 
shear  modulus,  respectively,  then  for  0  <  v  <  c*  the  stress  field  is  the 
same  as  for  elastic  material  (that  is,  it  is  independent  of  crack  speed 
and  material  properties),  whereas  for  c*  <  v  <  c  the  stress  field 
depends  upon  both  v  and  material  properties.  An  important  result  from 
this  study  was  that  the  dependence  upon  v  and  material  properties  has  a 
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simple  universal  form  from  which  qualitative  and  quantitative 
information  can  be  easily  obtained. 

In  this  paper,  the  corresponding  problem  for  a  viscoelastic  layer 
of  finite  thickness  is  considered.  The  finite  layer  problem  is  of 
interest  because  the  model  has  more  relevance  to  engineering 
applications  and  significantly  more  mathematical  complexity  than  does 
the  infinite  layer  problem.  Of  additional  importance  is  the  fact  that 
the  Riemann-Hilbert  method  employed  in  [2]  goes  through  for  this  more 
complicated  problem.  Specifically,  the  principal  obstacle  encountered 
in  applying  the  Riemann-Hilbert  method  is  the  evaluation  of  a  certain 
rather  complx  appearing  combination  of  functions  and  integral 
transforms.  It  is  demonstrated  here  that  a  simplification  similar  to 
that  effected  in  [2]  occurs  also  for  the  finite  layer  problem  and  that 
the  stress  intensity  factor  has  a  simple  universal  dependence  on  v, 
material  properties  and  layer  thickness,  even  for  general  shear  modulus. 

It  should  also  be  remarked  that  the  approach  adopted  here  has 
proved  to  be  convenient  for  the  calculation  of  the  angular  dependence  of 
the  near  crack  tip  stress  field  for  both  the  finite  and  infinite  layer 
problems.  These  calculations  are  the  subject  of  a  forthcoming  paper  as 
is  analysis  of  the  substantially  more  complicated  opening  mode  plane 
strain  problem.  This  section  concludes  with  the  formulation  of  the 
appropriate  boundary  value  problem  and  the  derivation  of  the 
corresponding  Riemann-Hilbert  problem. 

The  governing  field  equations  for  the  motion  of  the  linearly 
viscoelastic  layer  are 

°ij,j  =  P  ui,tt  '  —  <  *1  <  *  '  -h  i  x2  <  h  , 


°ij  =  2v*de^j  +  5£jA*dejcjt 

where  a^j,  e^j  and  u^  denote  the  stress,  strain  and  displacement  fields, 
respectively.  The  summation  convention  is  in  effect,  f,^  denotes 
partial  differentiation  of  the  function  f  and  p*de  denotes  the  Riemann- 
Stieltjes  convolution 

U*de  =  J  y(t-T)de(T)  . 

— oo 

Since  the  body  is  assumed  to  be  in  a  state  of  anti-plane  strain 
deformation,  the  only  equation  of  motion  not  identically  satisfied  is 

p*dAu3  *  p  u3>tt 

where  Au3  denotes  the  2-dimensional  Laplacian,  A  =  92/9x|  +  92/9 x|.  A 
semi-infinite  crack  is  assumed  to  be  propagating  to  the  right  with 
constant  speed  v  along  the  x^-axis.  The  crack  is  subjected  to  a 
traveling  distribution  of  applied  tractions,  o23(x1,0,t)  «  f(xj-vt)  for 
X1  I  vt,  while  on  the  upper  and  lower  surfaces  of  the  layer  two  possible 
boundary  conditions  are  considered,  fixed  grip  and  traction  free,  i.e. 

I.  u3(x^,±h,t)  =  0  ,  -®  <  x^  <  ® 

II.  °23 (xi,±h,t)  *  0  ,  -»  <  Xj^  <  ®  . 

Therefore,  adoption  of  the  Galilean  variables  x  =  x-^  -  vt,  y  =  x2  yields 
the  boundary  value  problem 

p*dAu3  *pv2u3^xx  |y|  <  h,  |x|  <  »  1.1 


o23(x,0) 


jy  (W*du3)  *  f(x) 


9 


x  <  0 
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I.  u3(x,th)  -  0  ,  | x |  < 


••  °23 (x»±h)  =  0  ,  | x |  < 


The  next  series  of  steps  are  similar  to  those  in  [2]  and  only  will  be 
summarized  here. 

Use  will  be  made  of  the  Fourier  transform  defined  by 
f(P#y)  =  j  e*pxf (x,y)dx  =  F+(p,y)  +  F_(p,y) 

— oo 

here 

F+  =  f  eipxf(x,y)dx  , 

•'O 

F_  =  /  e*pxf(x,y)dx  . 

hLoo 


Transforming  1.1,  solving  the  resulting  ordinary  differential  equation 
and  applying  the  boundary  conditions  1.3  yields 


u3(p,y)  =  A].<P) 


sinh(y(p)  (h-y)) 
cosh(Y(p)  (h-y)) 


with  Y ^  =  p*  +  ipv  /p(-vp).  As  will  be  discussed  later,  it  is 
convenient  to  choose  a  square  root  of  Y2  with  positive  real  part. 
Applying  the  boundary  conditions  1.2  produces  the  relation  on  y  =  0 


A  ,  A  A 

°23  +  °23  “  Wf-vpjY 


l(p)u3,l| 


coth  (hY(p) )  , 
tanh  (hY(p) )  , 


where  <*23  +  a23  denote  the  restrictions  of  o23  to  the  positive  and 
negative  x-axis,  repsectively. 

Equation  1.4  may  be  viewed  as  the  Riemann-Hilbert  problem 
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* 


F+(p)  =  G(p)F-(p)  +  g(p)  1.5 

where 

F+(p)  =  a+3 
F"(p)  =  U3#1 

A  A 

g(p)  =  -023  =  ~f 

*  Icoth  (hY (p) )  ,  I. 

G(p)  =  vU(-vp)Y(p)  -{ 

[tanh  (hY (p) )  ,  II.  . 

It  is  easily  demonstrated  a  posteriori  that  023(x,0)  and  U3^(xr0)  are 
such  that  F  (z)  =  o^3  and  F~(z)  =  uj  j  define  functions  analytic  for 
Im(z)  <  0,  respectively,  and  are  such  that  the  limits 

F±(p)  -  lim  F±(pfiq) 
q  +  01 

exist  and  satisfy  1.5.  In  the  next  section  the  Riemann-Hilbert  problem 
1.5  is  solved  and  the  stress  intensity  factor  calculated. 

2.  Problem  Solution 

In  order  to  solve  1.5  it  is  necessary  to  determine  the  mapping 
properties  of  the  coefficient  G(p).  (See  [1],  for  example,  for  a 
detailed  discussion  of  the  theory  of  Riemann-Hilbert  boundary  problems.) 
As  in  [2],  it  is  convenient  to  rewrite  G(p)  as 


where 


G(p)  *  sgn(p)G1(p)G2(p) 


Gx(p)  *  -iyt-vpjY^p) 


g2(p) 


coth  (hlpkj^p))  ,  I. 

tanh  (hJpjYjJp))  ,  II. 


2.1 


2.2 


2.4 


e“lvPt;du(t)  . 

o 

For  what  follows  it  is  sufficient  (but  clearly  not  necessary)  to 
assume,  as  in  [2],  that  the  shear  modulus,  u(t),  is  positive, 
continuously  differentiable,  non-increasing,  convex  and  such  that  y(°°)  = 
lim  y(t)  >  0.  These  assumptions  are  still  sufficiently  general  to 

t  “►<» 

include  all  physically  reasonable  examples.  Moreover,  an  easy 
adaptation  of  the  analysis  presented  here  permits  use  of  the  pure  power- 
law  model,  y(t)  *  uc(t/tc)-ot  ,  which  provides  an  effective 
characterization  of  marry  real  materials. 


For  convenience  the  following  observations  from  [2]  are  recorded 


here: 


(i)  y(o)  *  y(»)  <  Re(U(-vp))  <  u(o)  *  £(->) 


(ii)  lm(  jl(-vp))  =  -Im(y  (vp) ) 


>0,  p  >  0 

(iii)  arg(y(-vp) )  < 

[<0,  p  <  o 

(iv)  Im(Y^(-p))  =  lm(Y^(p)) 

Re(Y?(-p>)  =  Re(Y^  (p) ) 

(v)  Im(Yi<P))  >0,  0  <  p  <  « 

(vi)  1  -  (v/c*)2  *  yJ(0)  <  Re(Y2(p))  1  =  1  -  (v/c)2 

where  c*  *  (y («*>) / p)*/2  and  c  =  (w(0)/P)^/2  are  the  elastic  shear  wave 
speeds  corresponding  to  the  value  of  P(t)  for  infinite  and  zero  time. 


In  order  to  take  the  square  root  of  yf(p)  it  is  necessary  to 


distinguish  two  cases 

1.  O  <  v  <  c* 

2.  c*  <  v  <  c  . 

Taking  the  branch  cut  for  Y^(p)  to  be  the  negative  real  axis  assures 
that  Y-jJp)  has  positive  real  part.  Hence  for  case  1.,  y^  (p)  is  Holder 
continuous  for  all  real  p  and 


(vii)  Iin(Y1(p))  =  -Im(Y1(-p)) 

Re(Y1(p))  =  ReUjJ-p)) 

(viii)  ImfY^p))  >  0  ,  0  <  p  <  - 

(ix)  (1  -  (v/c*)2)1/2  =  y^O)  <  Y1(«.)  =  (1  -  (v/c)2)1/2  , 

whereas  for  case  2.  (vii)  and  (viii)  still  hold  but  Yx(p)  is 
discontinuous  for  p  =  0.  In  particular, 

(x)  Y^i®)  =  (1  -  (v/c)2)1/2 

yx(0 ±)  =  -  i ( (v/c*) 2  -  1)1/2  . 


The  image  in  the  complex  plane  of  the  real  p-axis  under  the 
transformation  Y^(p)  is  illustrated  in  Fig.  1  of  [2]  for  both  cases  1. 
and  2. 

Following  the  solution  method  employed  in  [2],  we  first  consider 
the  problem  of  finding  functions  X(z)  analytic  for  Im(z)  <  0, 

respectively,  and  which  satisfy  the  homogeneous  boundary  relation 


X  (p)  =  G(p)  X  (p)  . 


2.5 


Auxilliary  functions  X^z) ,  i  *  1,2  are  defined  by 


IKK' 


a‘V 


•  • « •  %  ■  *  *  %  *  • » •  • » * .  •  i 


Xj_(p)  =  Gt  (p)  XT(-p)  . 

Then  clearly, 

Xl z)  =  uT(z)  Xj(z)  X^z) 


2.6 


where  w  (z)  denote  branches  of  z^/2  whose  branch  cuts  are  the  negative 
and  positive  imaginary  axes,  respectively.  (See  [2].) 
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The  functions  X^(z)  were  constructed  in  [2] ,  but  here  only  the 
boundary  limit  x£(p)  is  needed.  It  was  shown  in  [2]  that 


(0 


+ 


(p)  x^p)  = 


(-ip) 1/2  Ig-l  («)  | 1/2  , 

(%-  ip) 1/2  |GX  (oo)  | 1/2  , 


where  cJq  is  the  unique  positive  solution  to 

vqo  r°°vi*(t)e“<3ovtdt  =  (v/c)2  . 
o 

In  2.7  u*(t)  denotes  a  normalized  modulus  given  by 


case  1. 
case  2. 


2.7 


M*(t)  =  y(t)/u(0)  . 


The  central  focus  of  this  paper  is  the  construction  of  the 
asymptotic  form  of  °23(x,0)  for  x  near  zero.  To  that  end,  much  of  the 
derivation  presented  in  [2]  is  valid  here  also.  In  particular,  it  can 
be  concluded  that 


F-(z)  =  X-tz)^ 


k( 


g(x)/X+(T)^ 


Jr  oo  />  oo 

e'ixPx+(P)dP5k  /  9(T)/X+(^)^ 

00  *^-oo 


and  for  x  near  zero,  the  dominant  term  in  the  asymptotic  expansion  of 
O23 (x)  is  given  by 


o23 (x)  a  x 


jz  I  g(T)/S(T)dx  =  k  x~l/z  . 


The  coefficient  K  is  the  stress  intensity  factor.  It  should  be  remarked 
that  2.8  is  valid  for  the  finite  layer  problem  because  G2(®)  =  1. 

In  [2]  a  simple  closed  form  expression  for  K  was  constructed 
through  the  device  of  introducing  the  function  h(x)  for  which 


'Z  (  T)  =  / 

“'-a 


e  ~1TXh(x)dx 


and  from  which  it  follows  that 


/CO 

°23 (x)h(x)dx  . 


The  principal  result  of  [2]  was  the  proof  that 


h(x)  =  j  G-^  (°°)  |  (-x)  |  x  | “1/2 


case  1. 


ex%  ,  case  2. 


where  H(x)  is  the  Heaviside  step  function. 


H(x)  = 


1,  x  >  0 

0,  x  <  0  . 


For  a  layer  of  finite  thickness,  lines  2.8-2.10  are  still  valid, 
.  .  _+ 

but  with  X  (t)  given  by  2.6.  Consequently,  the  simple  form  2.11  for 
h(x)  does  not  hold;  the  contribution  from  X^(x)  roust  be  incorporated 
into  the  solution.  It  is  convenient  to  remark  here  that  since  G2(p)  for 

-f 

problems  I  and  II  are  just  reciprocals,  the  same  is  true  for  ^(z). 
Therefore,  it  is  necessary  to  consider  only  one  of  the  cases,  say  I. 

Evaluating  X^(p)  is  the  central  problem.  Clearly,  from  the  general 
theory  of  Riemann-Hilbert  boundary  value  problems  ([1]),  it  follows  that 
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X*(z)  =  er2(z) 


’  2vr  / 


log  (G2  (t  ) )  . 


The  integral  in  2.12  can  be  evaluated  by  a  technique  similar  to  that 
employed  in  the  analysis  in  [2]  of  the  case  c*  <  v  <  c.  Moreover,  only 
the  boundary  limit  X2(p)  is  needed  in  the  computation  of  the  stress 
field.  As  in  [2],  it  is  necessary  to  consider  the  two  cases  1.  and  2. 


separately. 


Case  1. 


0  <  v  <  c* 


The  function  log  (G2(x)),  t  >  0  has  a  natural  extension,  log 

(G2(z)),  to  a  function  analytic  in  the  fourth  quadrant,  i.e.  for 

z  =  x  -  iq  with  x,q  >  0.  Computing  the  boundary  limit  lim  log 

X  +  0+ 

(G2  (x-iq) )  yields 

lhn^  log  (G2  (x  -iq) )  =  log  (coth  (-iqhYj^  (-iq) ) ) 

=  log  (i  cotdiqY-j^-iq)))  .  2.13 

Since  qY^(-iq)  is  an  increasing  function  of  q,  there  exists  a  sequence 
of  numbers,  (an}”=0,  such  that 

ha2nYl^“ia2n)  =  niT  '  n  =  0,1,... 

2.14 

ha2n+lYl(_ia2n+l)  =  (n+1/2)"  >  n  =  0,1...  . 


In  light  of  2.3,  2.14  can  be  rewritten 


va2n  f  y*(t)exp[-vta2n]dt  =  (v/c)2/[l-(mr/ha2n)2] 


va2n+l  f  P*(t)exp[-vta2n+1]dt  =  (v/c)2/[l-  (n+l/2)Tr/ha2n+1) 2] 


Consequently, 


arg  (i  cot  (hqY^-iq)))  =  */2  for  a2n  <  q  <  a2n+1  , 


while 

arg  (i  cot  (hqY-jJ-iq)))  =  -ir/2  for  a2n+1  <  q  <  a2(n+1)  . 

If  t  <  0  ,  the  function  log  (G2(t))  has  the  natural  analytic 
extension  log  (G2(z)),  z  =  t-iq,  q  >  0  given  by  log  (G2(z))  =  log 
(coth(-hz  Y^(z))).  Letting  z  approach  the  negative  imaginary  axis 
results  in 


THIo-  log  (G2  ( T-iq) )  «  log  (coth  (ihqY-^-iq)) ) 

=  log  (-i  cot  (hqY-^-iq)))  .  2.16 

The  function  r2(z)  may  be  computed  by  replacing  the  integral  in 
2.12  over  the  real  axis  with  the  appropriate  integrals  along  the  lines 
x"ia2n+l»  <  t  <  «  and  -iq,  0  <  q  <  a2n+1  .  Specifically,  for 


lm(z)  >  0, 


F2(z)  2rriX=c 


log  (G2(t))|I-2 


•1  f  r(0+)- ia^+l  f(0+)-io  r°>  , 

1  1  '  “  /  t  lcgIG  (,))i  2.17 

J(  0+) -ia^„ . i  JO  T  z 


0O-ia2n+l 


'(0+)-ia2n+1  JO 


1  f°  r(0-)-ia  2n+i  r~°°  ~  ia2n+l  ,  .  ..  dT  ^ 

orr  /  +  /  +  /  log  (g2  (t)  )  —x  2.18 

27rl  J-co  J(0-)-io  J(0-)-ia2n+l  2 


i  f  r(0+)-ia2n+i  r(0-)-io  , 

i  /  +/  log(G_(T))f^ 

2lT1  J(0+)-io  J(  0-)-iao_1,  2 


2.20 


1  f 00  ~  ia2n+l 
2ni  /  . 

-  ia2n+l 


log(G2(,))^ 


Lines  2.17  and  2.18  vanish  by  Cauchy's  theorem.  The  integral  in  line 
2.20  tends  to  zero  as  n  -*■  0  uniformly  in  z  for  Im(z)  >  0  and  therefore 
will  be  denoted  by  e(n). 

From  the  above  observations  and  after  an  obvious  change  of 
variables  in  2.19  it  is  now  clear  that 


r2(z)  =  e(n)  +JJJ 


n  Jo 


2n+1[log(icot(thY1(-it) ) )-log(-icot(thY1(-it) ) )]^§t 


1Z 


=  e(n) 


arg  (i<x>t  (thY-^  (-it) ) ) 


dt 

t-iz 


=  e(n)  + i-log(%^)-^[log(a: 

“1E=1 


In  2.21,  the  limit  Im(z)  -*•  0+  is  easily  taken  from  which  it  then  follows 
that 


-  (ip)"1/2  In(p)ee(n> 


2.22 
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with 

2n(P>  =  II  [  il^ip  ]  <a2n+l  "  iP) 1/2  •  2.23 

j=l 

Case  2.  c*  <  v  <  c 

In  contrast  to  case  l.#  for  case  2.  Y^(x)  has  a  natural  extension, 
Y^(z),  that  is  analytic  for  Im(z)  <  0  except  for  a  branch  cut,  -iq, 
0  <  <3  <  %  with  qQ  given  by  2.7.  The  same  is  then  true  for  log  (G2(x)) 
and  in  particular,  log  (G2(z))  is  analytic  for  Im(z)  <  -qQ.  After  this 
minor  modification  is  taken  into  account,  the  argument  utilized  in  case 
1.  applies  also  for  case  2.  and  yields  the  result 

X2(P)  =  (qo-iP)"1/2^(P)ee(n)  2.24 

where  ^(p)  is  still  given  by  2.23  and  the  an  by  2.14.  It  should  be 
observed  that  for  case  2,  the  an  are  such  that  qQ  <  a3  <  ...  whereas  for 
case  1,  0  <  a^  <  ...  .  Moreover,  the  an  for  case  2  are  larger  than  the 

corresponding  for  case  1. 

Calculating  a23(x,0)  requires  the  function  X.+  (p).  From  2.6,  2.8- 
2.10,  2.22-2.24  it  is  easily  seen  that  for  both  case  1  and  case  2 


X+(p)  = 


n 


lim 

n-+°° 


|g1( 


•>  i 1/2  TT  ( 

j=l 


aT_l  -1P 

%V’(a 


2n+l 


-  ip) 


1/2 


=  lim  X  +  (p) 
n->-«>  11 


2.25 


and  that 


K  *  lim  K_ 
n->°°  11 


,  =  -  ( Gx  (°°) 


1 1/2  £  f  ’2 


23(x,0)hn(x)dx 


with 


hn(x)  =  J^f  elxp/X„(p)dp. 


A  convenient  scheme  for  inductively  approximating  the  stress  intensity 
factor  K  can  be  based  upon  the  observation  that 


2.27 


where 


n  •  -  in  b  i  b 

J]  (-23,,^,  )=  x  +  nl  +  _  +  - 2" 

j=l  a2 j-1  -  lp  <arip>  (a2n-l"ip) 


bnk  =  [  J1  (a2-i  “  a2k-l)1/[  H  (a2j-l  _  a2k-l> J 

.  /  a2n  ~  a2k-l  N 

-  b(n-l)k'  a2n-i-a2k-1  ' 


In  light  of  2.25-2.27  and  after  some  routine  integrations  it  can  be 
shown  that  for  x  <  0, 


h^x)  *  |x|~1/2kn(x) 


kn(x)  =  ea2n-lx  -  2x 


l;  • 


exp[x(a2n+1t2  +  a2j-i (l~t2) ) Idt 
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